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Abstract 

We investigate one- and multi-dimensional stability of noncharacteristic boundary 
layers in the limit approaching a standing planar shock wave U{xi), xi > 0, obtaining 
necessary conditions of (i) weak stability of the limiting shock, (ii) weak stability of 
the constant layer u = U- := lim.z^-ooU{z), and (iii) nonnegativity of a modified 
Lopatinski determinant similar to that of the inviscid shock case. For Lax 1-shocks, we 
obtain equally simple sufficient conditions; for p-shocks, p > 1, the situation appears to 
be more complicated. Using these results, we determine stability of certain isentropic 
and full gas dynamical boundary-layers, generalizing earlier work of Serre-Zumbrun and 
Costanzino-Humphreys-Nguyen-Zumbrun. 
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Consider a boundary layer, or stationary solution, 

(1.1) U = U{xi), lim U{z) = U+, U{0) = Uq 

Z—t + OO 

of a hyperbolic-parabolic system of conservation laws on the quarter-space 

(1.2) Ut + Y,Fj{U)^^ = Y,{Bjk{U)Ux,).,, xeRi = {xi> 0}, t > 0, 

j jk 

with initial data U{x,0) = Uo{x) and boundary conditions as specified in ()1.6p below, that 
is noncharacteristic both in the hyperbolic sense 

det dF^{U+) / 

and with respect to the original (partially) parabolic problem as described in (H1)-(H3) 
below. Such layers occur, for instance, in gas- and magnetohydrodynamics with infiow or 
outfiow boundary conditions, for example in fiow around an airfoil with micro-suction or 
blowing; see |(;;MWZ51 IYZI iNZTj lNZ2] for further discussion. 
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As for any gas-dynamical flow, an important question is stability of these solutions 
under perturbation of the initial or boundary data. This question has been investigated 
in [SEl [MZn IGMWZ51 IGMWZ61 lYZl INZTI [NZ2] for arbitrary-amplitude boundary-layers 
using Evans function techniques, with the result that under quite general circumstances 
(see model assumptions below) linearized and nonlinear stability reduce to a generalized 
spectral stability condition phrased in terms of the Evans function, a Wronskian associated 
with the family of eigenvalue ODE obtained by Fourier transform in the transverse directions 
X := (x2, • • • , Xci)- See also the small-amplitude results of \GG\ \R3\ \MN\ iKNZl IKK] obtained 
by energy methods. 

The Evans function is readily evaluable numerically; see, e.g., [CHNZl [HLyZl HLyZ2| . 



As pointed out in [SZl ICHNZ[ IGMWZ5] . it is also evaluable analytically in certain inter- 
esting asymptotic limits. For example, it was shown in [GMWZSj that the Evans function 
converges in the small-amplitude limit as IJ approaches the constant layer U = [/+ to the 
Evans function of the constant layer, uniformly on compact sets of frequencies > 0, and, 
as a consequence, stability of small- amplitude layers is determined by stability of the limit- 
ing constant layer. This result was used in turn to show that noncharacteristic boundary 
layers of general symmetric-dissipative systems (defined below) are spectrally stable in the 
small- amplitude limit. 

A different asymptotic limit considered for special cases in \SZ\ ICHNZj is the standing 
shock limit X — >■ -l-oo in the case 

(1.3) U^{x) = U{xi - X), lim U{z) = U± 

2;— >±oo 

that U is the restriction to xi > of a standing shock solution U{- — X). It is natural 
to guess that there might be some relation between boundary-layer stability in this limit 
and stability of the limiting shock wave, and indeed it was shown in [CHNZJ for the case of 
isentropic ideal gas dynamics that the boundary layer Evans function, suitably normalized, 
converges in the standing-shock limit to the Evans function of the limiting shock wave on 
compact subsets of frequencies 3?A > 0, in complete analogy with the small-amplitude limit. 
On the other hand, it was shown in [SZj for the case of full (nonisentropic) ideal gas dynamics 
that boundary layers can in some parameter regimes be unstable in the standing-shock limit 
despite stability (see |HLyZl| ) of the limiting shock wave. 

In the present paper, we revisit the standing-shock limit in the general case, obtaining 
a result subsuming and illuminating these previous ones. Moreover, we carry out our 
investigations in multi-dimensions, whereas the analyses of [SZl ICHNZj were specific to 
the one-dimensional case. Specifically, we show that the boundary-layer Evans function, 
suitably normalized, converges in the standing-shock limit, uniformly on compact subsets 
of frequencies rj € M*^"^, 3?A > 0, to the product of the Evans function of the limiting shock 
wave and the Evans function of the constant layer U = U- at the left endstate of the shock. 
For symmetric-dissipative systems, this implies by stability of constant layers that the Evans 
function can be further renormalized so as to converge simply to the Evans function of the 
limiting shock, similarly as was shown in [CHNZ] in the special one-dimensional isentropic 
ideal gas case. 
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A consequence is that stability of both the hmiting shock and the constant layer U = C/_ 
are necessary conditions for stability of boundary layers in the standing-shock limit. On the 
other hand, these are not sufficient, as even stable shock waves have an (one-dimensional) 
eigenvalue at r/ = 0, A = due to translation-invariance, whereas stable boundary layers 
do not. A further necessary condition, therefore, is nonnegativity of the stability index 
(defined below) counting parity of the number of (one-dimensional) unstable roots rj = 0, 

> 0, indicating that the zero eigenvalue of the limiting shock does not perturb into the 
positive half-plane. Negativity of the stability index in the standing-shock limit, indicating 
an odd number of unstable eigenvalues, was what was shown in |SZj in order to obtain 
one-dimensional instability. 

Here, we develop these ideas substantially further, determining for Lax 1-shocks a simple 
and general stability determinant A(^, A,?]) extending the Lopatinski determinant A(^,A) 
of the inviscid shock case, G M, 77 G M, A € C, with 5RA, r] > 0, for which nonvanishing on 
the strictly positive half-space 5RA, r/ > is necessary and nonvanishing on the nonnegative 
half-space KA, rj >0 together with stability of the limiting shock U and the constant layer 
U = U- is sufficient for stability in the standing-shock limit, in one- and multi-dimensions. 
We then use this condition to investigate stability in various interesting situations. For p- 
shocks, p > 1, the situation is considerably more tricky, apparently involving a complicated 
double limit. 

1.1 Equations and assumptions 

Consider a family U^{x) of boundary-layers (11. 3p of (jl.2p consisting of translations of a 
standing shock solution U. Following [GMWZSl IGMWZ6] . we assume that the conservation 
law (jl.2p can be rewritten in nonconservative form, after an invertible change of variables 
U — > W, as a quasilinear hyperbolic-parabolic system 

d d 

(1.4) Ao{W)ut + Y,MW)djiW) - e dj{B,k{W)dkW) = 0, 

j=l j,k=l 

Aq invertible, with block structure 

(1.5) MW) = (I ;,,) , BMW) = (° , 

a corresponding splitting W = {w^,w'^) € R""*" x W, and decoupled boundary conditions 
(1-6) < '^2iw'^)\x&an = 92it,x), 

[ '^3{w,d:ri'w'^,ds:w'^)l^(,dQ = 0, 

where Ts{w,dx^w'^,dxw'^) = Kidx^w"^ + Y^'^j=2^ji'^)^XjW^ , Ki = constant, dimTi = n — r 
in the inflow case and in the outflow case (defined in (H2) just below), and 

dim T2 + dim T3 = r. 
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We make the following technical hypotheses following \Z1\ \Z3\ IGMWZS] . 
(HO) F^,B^^,A^,A^,Bj^,W{-),g{-,-) G C^ s > 2. 

(HI) The eigenvalues of ^j(^o^)~"^^]"'^Cj ^-re real and semisimple for all ^ 7^ in M'^. 

(H2) The eigenvalues of (^Q^)~^yl|^ are either strictly positive or strictly negative, that 
is, either a{Al^)-^Al^ > Oi > (inflow case) or (T(iji)-^i}i < -6*1 < (outflow case). 

(H3) ^(T^jkbi'^^j^k > 0\C\^ > for ah ^ G M"\{0}. 

(H4) The eigenvalues of dF^S^j are real, semisimple, and have constant multiplicity 
with respect to C&R'^,^^ 0. 

(H5) The eigenvalues of dF^{U±) are nonzero. 

(H6) iijdF^{U±) - ^■6^^'^(f^±)) < "^TW ^ ^ ^ > °- 

Definition 1.1. The system (jl.2p . (jl.4p zs symmetric dissipative at U± if in a neighborhood 
ofU± there exists a real matrix S{U) depending smoothly on U such that for such that for all 
^ G M'^\{0} S{U)Aq{U) is symmetric definite positive and block-diagonal, S{U) Ylj 
is symmetric, and ^S{U) ^ Bjk{U)^j^k nonnegative with kernel of dimension n — r. 

Alternative Hypothesis H4'. For systems that are symmetric dissipative at U±, we 
may relax (H4) to: 

(H4') About each ^ G M'^\{0}, the eigenvalues of dFj^S,j (necessarily real and semisim- 
ple, by symmetrizability) are either of constant multiplicity or else are totally nonglancing 
in the sense of jGMWZe], Definition 4.3. 

Definition 1.2. The system (|1.2p . (II. 4[) is symmetric hyperbolic-parabolic if there exists 
a real matrix S{U) depending smoothly on U such that for all ^ G M'^\{0} the matrix 
S{U)Aq{U) is symmetric positive definite and block-diagonal, {S{U)'^j Aj{U)S,jy^ is sym- 
metric, and the symmetric matrix ^S{U)'^ Bjk(U)£,jS,k is nonnegative with kernel of di- 
mension n — r. 

Examples 1.3. Hypotheses (H1)-(H6) are satisfied for standing shocks of the compressible 
Navier-Stokes equations with van der Waal equation of state, yielding boundary layers for 
which the normal velocity of the fluid is nonvanishing at Uq. This corresponds to the 
situation of a porous boundary through which fluid is pumped in or out, in contrast to 
the characteristic, no-flux boundary conditions encountered at a solid material interface for 
which normal velocity is set to zero. See |YZ1 IGMWZSl INZll INZ2| for further discussion of 
this situation and applications to aerodynamics. 

Hypotheses (H1)-(H6) with (H4) replaced by (H4') are satisfled for extreme (i.e., 1- or 
n-family) standing Lax shocks of the viscous MHD equations with van der Waal equation of 
state, with similar velocity restrictions on the plasma at Uq, but fail for intermediate shocks. 
Hypotheses (H1)-(H6) are generically satisfied for viscous MHD in dimension one, but fail 
always for viscous MHD in dimensions greater than or equal to two; see |MZ21 IGMWZSj 
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IGMWZ6] for further discussion. Both gas dynamics and MHD equations with van der Waals 
equation of state are symmetric hyperbohc-paraboUc systems that are symmetric dissipative 
at U± for standing shocks connecting thermodynamicahy stable endstates jZ3l IGMWZ4j . 

Finahy, regarding the standing shock U, ordering the eigenvalues of dF(U±) as 

af <af-- - <a^, 

we assume: 

(H7) Profile C/ is a transversal viscous Lax p-shock, i.e., 

(1.7) OpLi < < Op , a+ < < a+_^ 

and [/ is a transversal connection of the standing wave ODE with boundary conditions U± 
(see |MaZ3| for a detailed discussion of the standing wave ODE) . 

The eigenvalues correspond to characteristic speeds at U± of the associated one- 
dimensional inviscid system Ut + Fi{U)xi = 0, with ap the principal characteristic speed 
associated with the shock. 

1.2 The Evans condition 

The linearized eigenvalue equations of (jl.2p . (jl.6p about U are 

(1.8) XU = LU := 5^(i?,feC/.J., - Y.^A,U)x^, 

j,k 3 
with homogeneous boundary conditions 

(1.9) r{w,diw\dx,w\,^=o = ^ 

expressed in linearized W-coordinates W := {dW / dU){U)U , where W and U denote per- 
turbations of W and U. 

Taking the Fourier transform in x := (x2, . . . , x^), we obtain a family of eigenvalue ODE 

LoU 

XU = L^U := {BiiUy - {AiUy -i AjCjU + i ^ Bji^jU' 

(1.10) m m 

with boundary conditions 

(1.11) T'{W,iiw^,dx,w%^=o = 0- 
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1.2.1 The boundary-layer Evans function 

A necessary condition for linearized stability is weak spectral stability, defined as nonexis- 
tence of unstable spectra 5RA > of the linearized operator L about the wave. As described 
in Section [2l this is equivalent to nonvanishing for all € M'^"^, 5RA > of the Evans 
function 

D{1\) 

a Wronskian associated with (jl.lOp with columns consisting of bases of the subspace of 
solutions decaying as xi — > +00 and the subspace of solutions satisfying the boundary 
condition (jl.lip . Under our hypotheses, the Evans function may be defined to be C°° away 
from the origin on > with continous limits (typically depending on direction) at (0, 0) 
along rays through the origin; see |GMWZ5l IGMWZ6] . 

Definition 1.4. We define strong spectral, or uniform Evans stability as 

(D) \D{lX)\>e[C)>d 

for (I, A) on bounded subsets C C G W^-^ , KA > 0} \ {0}. 

Remark 1.5. Under assumptions (H0)-(H6), uniform Evans stability implies linearized 
and nonlinear stability in both the long time and small viscosity limits of general nonchar- 
acteristic boundary layers (not necessarily associated with standing shocks) of symmetric 
hyperbolic-parabolic systems that are symmetric-dissipative at C/+ with Dirichlet boundary 
conditions dimTs = 0; see |GMWZ5[ iGMWZGi INZTI lNZ2l [N2] . For more general systems 
and boundary conditions, (D) augmented with a rescaled high-frequency condition has been 
shown in [GMWZSl IGMWZ6] to imply stability in the small viscosity limit. 

1.2.2 The shock Evans function 

Likewise, a necessary condition for linearized stability of the shock wave tl is weak spectral 
stability, defined as nonexistence of unstable spectra 3?A > of the linearized operator L 
about the wave, or nonvanishing for all ^ E M*^"^, JRA > of the shock Evans function 

a Wronskian associated with (jl.lOp with columns consisting of bases of the subspace of 
solutions decaying as xi — > -|-oo and the subspace of solutions decaying as xi — > -co. 
Under our hypotheses, the shock Evans function may be defined to be C°° away from the 
origin on KA > and C° at the origin, with first directional derivatives (typically depending 
on direction) at (0,0) along rays through the origin; see [Z3 l IGMWZ41 IGMWZ6] . 

Definition 1.6. Uniform Evans stability of a standing shock is defined as 

{V) \v{t\)\>e{c)>Q 

for (I, A) on bounded subsets C C {| G W^'^, KA > 0} \ {0}. 
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Remark 1.7. Under assumptions (H0)-(H7), uniform Evans stability (P) implies lin- 
earized and nonlinear stability in both the long time and small viscosity limits of standing 
shocks of symmetric hyperbolic-parabolic systems that are symmetric-dissipative at U±; 
see |GMWZ4[[Z3l lN2]. 

1.3 Main results 

I. 3.1 Convergence 

Denote by Dx{i, A) the Evans function associated with tj-^ and D^{^, A) the Evans function 
associated with the constant boundary-layer U = C/_ at the lefthand endstate C/_ of U. 
Then, our first main result is as follows. 

Theorem 1.8. Under assumptions (H0)-(H7), there exists a continuous nonvanishing 
function (3{^, A, X) such that 

(1.12) DxilX) ^ l3{iX,X)D^{i,X)V{lX) 

as X ^ oo, uniformly on compact subsets of € M"'^"'^, KA > 0}. 

Corollary 1.9. Under (H0)~(H7), weak spectral stability of both the constant boundary- 
layer U = U- and the standing shock U are necessary conditions for stability of U^ in the 
standing-shock limit X — > oo. 

Remarks 1.10. 1. Under the stronger definition of uniform Evans stability defined in 
[GMWZSl IGMWZ6j involving also a rescaled high-frequency condition, and assuming (HO)- 
(H7), uniform Evans stability of the constant boundary-layer U = f7_ and the standing 
shock U are also sufficient conditions for spectral stability of in the standing-shock 
limit X ^ oo for frequencies uniformly bounded away from the origin (^, A) = (0,0). For 
intermediate frequencies R^^ < |(^, A)| < R, this is an immediate consequence of Theorem 

II. 8i For high frequencies |(^, A)| > i?, i? > sufficiently large, it follows from the fact 
established in Section 3.2 [GMWZ5] that high-frequency stability is equivalent to stability 
of the constant layer U = Uq, and the fact that Uq^ — > C/_ as X ^ oo. That is, assuming 
stability of U and U = [/_, unstable frequencies, if they occur, must converge to the origin 
as X ^ oo, with no additional assumptions on the system ()1.2p . 

2. It is shown in Corollary 1.29 [GMWZ5] that constant boundary-layers of symmetric- 
dissipative systems are uniformly Evans stable (under the stronger definition of }GMWZ5t 
IGMWZ6] ) for Dirichlet boundary conditions. Thus, if (jl.2p is symmetric-dissipative at U- 
and the boundary conditions are Dirichlet-type, then the constant layer U = U- is stable, 
with the implications above. In the general case, stability or instability of the constant layer 
U = U- may be determined by a linear algebraic computation, since the constant-coefficient 
eigenvalue ODE is explicitly soluble for each frequencies (^, A); see [GMWZS] for further 
discussion. 
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1.3.2 One-dimensional stability 

Due to translational invariance, T>(0, 0) = 0, and so we cannot conclude nonvanishing of Dx 
near the origin as X — > oo from the convergence result (jl.l2p . Indeed, it is possible that the 
zero of V at the origin may perturb into the unstable half-plane 5ftA > for boundary layers 
with X large, yielding instability. In the one-dimensional setting, this may be dectected by 
the stability index 

r := sgn lim D{0,X)x lim sgnZ)(0, A), 

A^o+real A^+cx)real 

where D is chosen with a standard normalization guaranteeing that it is real for real A and 
^ = 0. The stability index is well-defined by the properties that D is continuous along 
rays at the origin and nonvanishing for real A sufficiently large; see IGZ\ \SZ\ IZ3| for further 
discussion. Negativity of F, by the Intermediate Value Theorem, implies existence of a real 
positive root I?(0, A*) = 0, hence one- dimensional instability. 

Observing that stability of the constant layer U = C/_ and the limiting shock U imply 
that sgnD_(A) and sgnX'(A) are constant for real A > 0, we thus obtain by convergence, 
(jl.l2p . that a necessary condition for stability in the standing shock limit X +oo, 
assuming stability of the constant layer, is nonnegativity of 

(1.13) f:= lim sgn lim tT^t^ ^ 1™ sgnP(0,A), 

x^+oo A^o+real -C-(0, A) A^o+real 

provided these limits exist. 

A standard shock stability result |GZl \ZS\ IZSj is that, with appropriate normalization, 

lim sgnP(0, A) = sgndx'D{0, 0) = sgn.6, 
(^1.14) A^o+real 

6 := sgndet(i?",i?+, [U]), 

where and are matrix blocks whose columns span the stable subspace of dFi{U^) 
and the unstable subspace of dFi[Uj^) and \U] := C/+ — U- denotes the jump across the 
shock. The determinant det(i?~,i?+, \U]) may be recognized as the Lopatinski determinant 
of one-dimensional inviscid theory, whose nonvanishing is equivalent to one-dimensional 
inviscid shock stability. 

Our second main result asserts that the first limit in (jl.lSp also exists, yielding a neces- 
sary stability condition of nonnegativity of a certain Lopatinski-like determinant 5 relative 
to the sign of 6. 

Theorem 1.11. Assuming (H0)-(H7) and stability of the constant-layer U = U- and the 
limiting shock U, 

lim sgn lim = sgno, 

x^+oo x^o+ real D- (0, A) 

,5 := sgndet(i?",fl+,F) 
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provided 5^0 and S := lim^^_oo exists, where R has the dimensions of R and 

V is a single column vector, hence 

(1.15) sgn56 = sgndet{R-,R+,[U])det{R-,R+,V) >0 

is necessary for one-dimensional stability in the standing- shock shock limit. 
For a Lax 1- shock with Dirichlet boundary conditions, 

(1.16) 5 = det(i?+, [U]) and 5 = det{R+,dFi{U^)S) 
so that (ll.lSp becomes 

(1.17) sgndet(i?+, [U])det{R+,dFi{U^)S) > 0. 

Remarks 1.12. 1. In the (characteristic) hmit C/- ^ C/+ as the amphtude of the back- 
ground shock U goes to zero, S* ~ [U] ~ r^, where ri is the eigenvector of Ai associated 
with the smallest eigenvalue ai, and dFi{U-)S ~ A'^r^ = aj~[C/], where > by the Lax 
shock conditions dTTl)- Thus, sgndet(i?+, [U]) dei{R+, dFi{U_)S) ~ sgndet(i?+, [C/])^ > 0, 
consistent with stability. 

2. Prom the boundary-layer ODE, dFi{U)U' = {Bn{U)U'y, giving Ai{xi)U' = BnU", 
where ' denotes S^i- Thus, = aB^^S, where a = limz^-oo{U' /U){z) is necessarily real 
and positive if the limit S := lim2^_oo(f/'/|C/'|)(2:) exists. It follows that we may replace 
(jl.l7p by the equivalent condition 

(1.18) det{R+,[U]))det{R+,Bn{U-)S) > 0, 

which may be recognized as the necessary condition derived by a rather different argument 
in Section 4.1 of [SZj. 

Our third result states that for Lax 1-shocks the necessary conditions we have derived 
are also essentially sufficient for one-dimensional stability in the standing shock limit. 

Theorem 1.13. For a Lax 1-shock with general boundary conditions (11.60 . assuming (HO)- 
(H7), positivity ofT together with stability of the constant layer U = and the limiting 
shock U is sufficient for stability in the standing-shock limit X -\-oo. 

Remarks 1.14. 1. It was shown in [HuZj that in the limit U- f7+, the background shock 
U is stable. For symmetric-dissipative systems with Dirichlet boundary conditions, the con- 
stant layer U = C/_ is stable. By Remark 11.121 1. therefore, for Lax 1-shocks of symmetric- 
dissipative systems, with Dirichlet boundary conditions, layers are one-dimensionally 
stable in the standing-shock limit X — > -|-oo for U-^- fixed and shock amplitude |{7+ — U-\ 
sufficiently small. 

2. By Remark 11.121 2. the necessary condition of [SZj . together with stability of the 
limiting shock U' and the constant layer U = U-, is also sufficient for one-dimensional 
stability. For ideal gas dynamics, the numerical study of |HLyZlj indicates one-dimensional 
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stability of arbitrary shock waves U for gas constant 7 within the physical range 1.2 < 7 < 3 
(the only values considered); likewise, U = U- is stable by symmetric-dissipativity of the 
compressible Navier-Stokes equations. Thus, for ideal gas dynamics with Dirichlet boundary 
conditions and 1.2 < 7 < 3, one- dimensional stability in the standing Lax 1-shock limit is 
completely decided by the simple algebraic condition (|1.18p of JSZf . 

For Lax p-shocks, p > 2, the situation is more complicated, apparently involving a 
tricky double limit. In particular, the conditions of Theorem 11.111 are only necessary and 
not sufficient for stability. 

1.3.3 Multi-dimensional stability 

We restrict now to the case of a Lax 1-shock, for simplicity taking pure Dirichlet boundary 
conditions, dimTa = 0. In multi-dimensions, uniform inviscid stability of a Lax 1-shock is 
defined as nonvanishing of the multi-dimensional Lopatinski determinant 

d 

(1.19) A(f , A) := det (7^+(|, A), X[U] + iCjlFjiU)]) , 

i=2 

on the nonnegative unit sphere := {[(C, A)| = 1, KA > 0}, where TZ^ is a matrix blocks 
whose columns form a basis for the unstable subspaces of 

d 

(1.20) A+il A) := (a/ + J2i^jdF,{U+)yFi{U+r' 

i=2 

and [/i(f/)] := h{Uj^) — h{U-) denotes jump in h across the shock. Define the related 
determinant 

d 

(1.21) A(e,A,r/) ■.= det(n+ilX),X[U] + ^iC,[Fj{U)]+rjdFi{U.)S.), 

i=2 

E M, where S := lim^^-oo j^(-z). Then, our fourth and fifth main results, giving neces- 
sary conditions and sufficient conditions for multi-dimensional stability analogous to those 
of Theorems 11.111 and 11.131 ^-re as follows. 

Theorem 1.15. For a Lax 1-shock and Dirichlet boundary conditions, assuming (HO)- 
(H7), a necessary condition for stability ofUx in the standing- shock limit X +00 is that 
A have no root that is simple with respect to A on the positive half-sphere 

:={|e,A,r?| = l, KA>0, ry>0}, 

in the sense that A = and d\A ^ 0. 



11 



Theorem 1.16. For a Lax 1-shock and Dirichlet boundary conditions, assuming (HO)- 
(H7), sufficient conditions for stability of Ux in the standing- shock limit X +00 are 
nonvanishing of A on the nonnegative half-sphere := {|^,A,r/| = 1, 3ftA > 0, r] > 0}, 
stability of the limiting shock U , and stability of the constant layer U = 

Remark 1.17. In the one-dimensional case ^ = 0, TZ± = R± = constant, so that d\A = 5. 
It is readily seen that Theorems I1.15l and ll.l6l reduce in this context to the restrictions of 
Theorems 11.111 and 11.131 to the case of a Lax 1-shock and Dirichlet boundary conditions, 
provided that the limiting shock U is one-dimensionally stable, so that 5^0 [ZH] . 

Remarks 1.18. 1. Multi-dimensional stability of shock waves in the limit — > [/+ has 
been established in |FS| for symmetric dissipative systems with strictly parabolic (Lapla- 
cian) viscosity Yli'^jkUxj)xk By the arguments of Remarks 11.121 1 and 11.1^ 1. therefore, 
boundary layers of such systems with Dirichlet boundary conditions are multi-dimensionally 
stable in the standing Lax 1-shock limit for U-^- fixed and — U- \ sufficiently small. The 
argument of [FS] appears likely to generalize to the general symmetric dissipative case (see 
also [PZ]), which would extend the boundary layer result also to the general case. 

2. For ideal gas dynamics, the numerical study of [IlLyZ2| indicates multi-dimensional 
stability of arbitrary shock waves for gas constants in the physical range 1.2 < 7 < 3. By 
the arguments of Remarks [LT2j 2 and ll.l4l 2. therefore, for ideal gas dynamics with Dirichlet 
boundary conditions and 1.2 < 7 < 3, multi- dimensional stability in the standing Lax 1- 
shock limit is completely decided by vanishing or nonvanishing of the extended Lopatinski 
determinant A defined in ()1.2ip .' a simple, linear-algebraic condition. 

1.3.4 Verification 

Evidently, nonvanishing of A on the nonnegative half-sphere in (^, A, rj) is equivalent to the 
condition that the image of 

(1.22) fiilX) := -A(e, A)/det (7^+(|, A), dFi(C/_)5_) 

over the nonnegative half-sphere in (^, A) avoid the nonnegative real axis, a natural gener- 
alization of the one-dimensional condition (jl.lTp . Note that rj is independent of the choice 
of 71+ and homogeneous degree one in (^, A). This leads us to the following condition 
convenient for numerical or analytic verification. 

Proposition 1.19. Assuming the one- dimensional stability condition ()1.17p . nonvanishing 
of A on the nonnegative half-sphere in (^,A, ry) is equivalent to the condition that the image 
of f}{\,iT) over r G R, with fj defined as in ()1.22p . does not intersect the nonnegative real 
axis. 

Proof. By homogeneity of A, nonvanishing on the half-sphere is equivalent to nonvanishing 
of (5(0, A, ry), or (|1.17p . and nonvanishing of (5(1, A, rj), or the condition that r/(l. A) avoid the 
nonnegative real axis for 5RA > 0. Recalling the standard fact that 7^+(l, A), hence 77(1, A), 
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may be chosen to be analytic in A for 5RA > and continuous at ?R.X = 0, we find by the 
argument principle applied to a sufficiently large semicircle about the origin, bounded to 
the left by the imaginary axis, and the fact that 77(1, A) by homogeneity /continuity does 
not intersect the nonnegative real axis for KA > and |A| sufficiently large (by the assumed 
one-dimensional stability) that the latter condition is equivalent to nonintersection with the 
nonnegative real axis for as A traverses the imaginary axis. □ 

1.4 Discusssion and open problems 

The results of Theorems 11.151 and 11.161 illuminate and greatly extend the earlier results of 
[SZ] and [CHNZ] in the one-dimensional case. In particular, we regard the derivation of 
necessary and sufficient conditions for multi-dimensional stability as a substantial advance. 
Though our necessary and our sufficient conditions are slightly different, the difference 
is sufficiently slight that it should not interfere in practice with classification of physical 
stability regions. 

We note that, besides its independent interest, the treatment of the standing shock 
limit, as pointed out in [CHNZ] . is important in truncating the computational domain for 
global stability analyses. 

On the other hand, we have restricted here mainly to the simplest case of a Lax 1-shock, 
which corresponds to the case of inflow boundary conditions. It would be very interesting 
to obtain corresponding conditions also in the case of outflow boundary conditions, for 
example, a Lax n-shock. Likewise, it would be very interesting to carry out computations 
analogous to those carried out for gas dynamics in Section [5] also for the equations of MHD. 

2 Construction of the Evans function 

We begin by reviewing the construction of the Evans function following \Z3\ IGMWZSj 
l(^MWZ6[[Al[NZ2] . 

2.1 Expression as a first-order system 

We first observe that matrix 



is full rank, by (H2)-(H3) together with block structure assumption (jl.4p . as can be seen 
most easily by working in M^-coordinates; see |Zll IZ31 [MaZSj . 

As a consequence, the Fourier-transformed eigenvalue equations (jl.lOp . (jl.lip may be 
written as a first-order system 




(2.1) 



Z' = Gxii,X,xi)Z, 
M{i)Z = 



13 



in the convenient coordinates 



(2.2) 



Z 




BuU' -AiU 



( Bflu + Blfv \ 

-Al^u-A\^v 
Bllu' + Bllv' - Afu - Afv 



\ 



J 



see |Z3] for further discussion. 

In the case of Dirichlet boundary conditions dimTs = 0, the boundary operator M is 
independent of ^, with 



(2.3) 

in the inflow case and 
(2.4) 



ker M = {Z : Zi = 0, Z2 = 0} 



kerM = {Z : Zi = 0} 



in the outflow case. In general, keri?(^) depends on ^ in a possibly complicated way. By 
the underlying relation Ux{xi) = U{xi — X), we have 



(2.5) 



Gx(|, A,a;i) 



, A,rci - X). 



2.2 Conjugation to constant-coefficients 



We next recall the following result established in [Z31 IGMWZ5] . a consequence of the 
conjugation lemma introduced in [MZlj and the fact proved in |MaZ31 [Z3j that U under 
hypotheses (H0)-(H7) converges exponentially to U± as xi — > ±00. See, e.g., |Z3] for details. 



Proposition 2.1. There exist matrix-valued functions 

T^{iX,xi) = T^{lX,xi-X), 
uniformly bounded with bounded inverse for xi ^ 0, locally analytic in (C,A), such that 
(2.6) |(T± - Id){xi)\ < Ce-^'l^il for xi ^ 0, 

9 > 0, and Z := T^X satisfies the constant- coefficient equation 



X' 



(I, X)X for xi^X 



G-^{^,X) := Gj^(^, A, ±00), whenever Z satisfies (12. ip . 
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2.3 Definition of the Evans function 

Finally, we recall the following standard result, established in increasing generality in \SZ\ 
IZTj Izi l IGMWZI IGMWZ6] . 

Proposition 2.2. Under (H0)-(H7), there exist matrices E^, E+, and E'^ whose columns 
form bases of the unstable subspace of G~ , the stable subspace of , and ker7V4(,^) and 
matrices ¥~ and whose columns form bases of the the stable subspace of G~ and the 
unstable subspace of , C°° on {^X > 0} \ {(0,0)} and continuously extendable along 
rays through the origin. Moreover, E"*" and F+ and E~ and F^ are uniformly transverse on 
compact subsets of {^X > 0} and 

(2.7) dim Span E^ = dim Span ¥P = {n -\- r) — dim Span E^. 

Remark 2.3. In the one-dimensional setting, the subspaces of Proposition 12.21 may be 
defined as globally analytic functions on $RA > —r], rj > 0, using a standard construction of 
Kato [lit]; see, e.g., [Z3l ICHNZ] . 

With these preparations, we define the shock Evans function precisely as the (n + r) x 
(n + r) Wronskian 

(2.8) V{1 A) := det (t-E", T+E+) U,=o 
and the boundary-layer Evans function as 

(2.9) DxilX) ■■= det (e", T+E+) U.,=o = det (e", r+E+) U,=_x. 
The Evans function for the constant-layer U = U- is given by 

D^{i,X) := det (e°,F-), 
since in this case E+ = F^, or equivalently by 

(2.10) D-{i,X) := det (^E-*E°) , 
where 

(E-,F-) = (E-,F-)-i* 

are dual bases to (E^, F^) with respect to the standard complex inner product, M* denoting 
adjoint, or conjugate transpose, of a matrix M. 

2.4 Behavior near zero 

For later use, we record the following refinement of Proposition 12.21 (also established in 
[pl EH ESI l(;;MWZ4llGMWZ6p . from which we may determine the behavior of V, as 
(|,A)^(0,0). 
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Proposition 2.4. Under (H0)~(H7), the limiting subspaces 

lim Span T+e'^^''^E+(pio, pXo) and lim Span e^'^^'E- (pio, pXo) 

/9 £ M, KAo > 0, are spanned by the direct sum of fast modes 

Zi = bii(t), (Z2,Z3) = 0, 
with (j) satisfying Bii<j)' — Aicf) = and decaying at +oo [resp. —oo] and slow modes 

Zi = *, (Z2,Z3) = 7^, 

with TZ spanning the unstable [resp. stable] subspace of [resp. A^, with 

d 

A±{lX) := (m + Y.iCjdFj{U±)yFi{U±r\ 
i=2 

expressible alternatively as Z = T^X with 

(2.11) Xi = bu{-Af)-^n, {X2,Xs) = n 

constant. Symmetric decompositions hold for T^e^^^^¥^ andT~e^ ^^F". 

3 Basic convergence result 

Proof of Theorem \1.8\ for A) bounded away from the origin. Viewing V and Dx as Wron- 
skians of solutions of the same ODE (|2.1|) , we may rewrite (|2.9|) using Abel's Theorem as 

(3.1) Dx(,i,X) ■■= ei'o"''TraceG«,A,.)d^^g^ f 5-^-"°E°, T+E+) ^,=0, 



where S^^^ denotes the solution operator of (j2.ip . 
Next, expand 

E° = r~(-x)r"(-x)"^E° 

(3.2) =T-E%,=^x + 0{e-''') 

= T-U^-E%,=_x + T~U^-EX,=_x + 0(e-'^), 

where 

(3.3) He- = E^E^* and Hp- = F^F"* 

denote the eigenprojections of onto subspaces E~ and F~, noting that 

(3.4) r"niE-E° = (r"E")(E"*E°). 
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From (13.41). we obtain 



(3.5) = det (S-^^^T-E-,T+K+j |^,=o det(E-*E°) 

= det (s-^'^^T-E-,T+E+^ Ui=o-D_. 

Next, expanding 
and noting that 

(3.6) e'^'^^'-y^E- = enE-'G-{^-s/)n]g_E- = E-e^"''^"'^"^^"^'), 
where Trace(E-*G"E") = Trace(E"*E"G") = TracenE-G", we find that 

det(5"^^°r-E-,r+E+)|^.,=o 

(3.7) ^gTrace(n^_G-)x^g^ (^-E-,^+E^ 

^ gTrace(nE_G-)Xp 



where, since E is the unstable subspace of G , 

gTrace(nj._G-)X 

is uniformly exponentially growing in X for 3?A > 0. Indeed, for (^, A) bounded from the 
origin 5ftA > 0, each column of S^^~^^T~E~ is uniformly exponentially growing in X at 
rate at least 

where fi*{6., A) is the smallest real part of the (positive real part) eigenvalues of G^ associated 
with the unstable subspace E~. 

By a similar argument, each column of 5~"''^^'^r~nip'-E'' is uniformly exponentially 
decaying in X, at rate e^*-^ , where ^* is the largest real part of the (negative real part) 
eigenvalues of G~ associated with the stable subspace F~. Collecting information, we thus 
have 

Dxil A) = m, A) [DxH \)V{1 A) + 0(e-^^) + 0(e(^*-'^*)^)) 
^ (i{l\X)Dx{lX)V{l A) 
as X — > +00, exponentially in X, where 

for (^, A) uniformly bounded away from the origin on KA > 0. □ 
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Remark 3.1. In general, the real parts of the eigenvalues of E~ and F~ can converge to 
zero as A) approach the origin, so that /x*, /x* ^ and the above convergence argument 
fails. In the special case of a Lax 1-shock, /i* 0, but /i* remains strictly negative [Z3], 
and so we obtain uniform convergence by this argument on all of > 0. Indeed, in 
the one-dimensional case, we obtain uniform convergence on 5RA > —rj, ry > 0. (The only 
obstruction in the multi-dimensional case is that D is not defined on this set [Z3j.) 

Combining Remarks 11.10^ 12.3^ and 13. H we obtain the following simple result reducing 
determination of one-dimensional stability to computation of the stability index. 

Lemma 3.2. Under (H0)-(H7) and stability of the constant layer U = C/_, for a stable Lax 
1-shock in one dimension, Dx has exactly one zero on 5RA > —r], rj > 0, for X sufficiently 
large, hence is stable if and only if its stability index is positive. 

Proof. By Remarks I1.10( under stability of C/ = U- , we may restrict attention to a compact 
set of frequencies, while by Remark 1 2. 3 1 we may take Dx and analytic on 5RA > —rj. As the 
uniform limit of analytic functions, we find that the number of zeros of VD- on 5RA > —rj 
is equal to the number of zeros of Dx for X sufficiently large. As the number of zeros of 
P is one and the number of zeros of D^ is zero, by stability, the number of zeros of Dx is 
one as asserted. 

As the stability index T counts the parity of the number of nonstable roots 5RA > (see 
\GZ\ IZ3| ). positivity of T corresponding to an even number of nonstable roots, and since 
r = corresponds to instability [NZlj . we thus obtain stability if and only if F > 0. □ 

4 Behavior near the origin 

Proof of Theorem \1.8\ for (^, A) (0, 0). The shock evans function D is continuous at the 
origin, with V{0, 0) = 0, and Z)_ is bounded on compact sets. Thus, to complete the proof 
of Theorem 11.81 it suffices to show that 

Dxil A)M, A)^0 

as (^, A) — > and X — > +oo, with /? defined as in (j3.9p . 

Noting that T^E^(0, 0) contains by continuity all exponentially decaying solutions of 
the one-dimensional eigenvalue equation, hence, in particular, U'{xi), we may without loss 
of generality assign the value U'{xi) to the first column of T~E~ at (0,0). Moreover, 
noting that the strongly unstable subspace of , defined as the part whose eigenvalues 
have strictly positive real part even at (0,0), perturbs analytically, we may restrict the first 
column to this subspace, ensuring that the first column of S~^^^T~E,~ is analytic up to 
the origin and moreover grows exponentially in X, at rate e^^, some 9 > 0, for |(^,A)| 
sufficiently small, hence contributions to Dx coming from the first columns of the second 
two terms in the last line of (j3.2p are exponentially small as X ^ and can be ignored, 
while the contributions in other columns are at least bounded. Likewise, we may arrange 
that the first column of be analytic up to the origin and equal to U' at (^, A) = (0, 0). 
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In place of (j3.8p . therefore, we obtain the weaker estimate 



l^x(e,A)|//5(C,A) < Dx(C,A)det T-(Er,0(l)),T+E^ 

\ / xi=0 

(4.1) +0(e-^^)+0(e(^-^*)^) 

^ 

as A) — > (0,0), X —5- +00, where T^E^^ and T^K'^ denote the first columns of T^E^ 
and T+E"*", since by our choice of normalization T^E^^ and T+E^ are continuous (indeed, 
analytic) at the origin and coincide for (^,A) = (0,0). This completes the proof of the 
theorem. □ 

Proof of Theorem Fix ^ = 0, so that -D(0, A) is continuous in A. Without loss of 

generality, take Z)_(0,0) = 1, with, moreover, 

(4.2) E°(0,0)e-'^"*'^"'^"-^ = r(-X)e-^"^(^E-(0,0) +F-(0,0)a) 
for some k x k matrix a, where k = dimE~. Thus, 

^-x^o^o ^ T'{0)(e'^'^T-{-X)-^E'^' 



(4.3) 

= T-(0)(e~(0,0) +E'(0,0)a)e-^ *® ^ ^. 



1x1=0, (|,A)=(0,0)' 



Following the proof of Theorem 11.81 we find that the quantity 

Dx{0, 0) /C'^o"^ TraceG(0,0,2;)d2gTraceng_G-X 

is given (exactly, with no exponentially decaying error) by 

det (^T" (E" +E"a, r+E+^ 
where TraceG(0, 0, z) and Tracen£;£;-G~(0, 0) are real, hence 

g/g"^ TraceG(0,0,2)d2 gTraceng_G~X 

are real and positive. 

Appealing to Proposition 12.41 we may arrange without loss of generality that the first 
ki columns of T~(0, 0, xi)e'^ (°''^)^'iE~(0, 0) consist of functions {bii(j)j,0), where j = 
1, . . . ,ki are solutions of 

(4.4) 5n0;. - ^10, = 

that are uniformly exponentially decaying as xi — oo, hence exponentially growing in 
forward direction, and the remaining k2 columns consist of functions (*, r~), j = 1, . . . , /c2, 
where r~ are constant eigenvectors of with negative eigenvalues aj, and similarly for 
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r+E+. Likewise, we may arrange that the columns of T F consist of li solutions of (j4.4p 
that are uniformly exponentially decaying in xi in forward direction and I2 functions {*,rj), 
j = 1, . . . ,l2, where rj are constant eigenvectors of with positive eigenvalues aj. 

Finally, we may choose the first column of T~'E~\xi=o as U'{0), noting that U'{0) lies 
also in T"'"E"'"|a;^=o- Combining these facts, we find that, up to an exponentially decaying 
error with respect to X, we may rewrite Z)x(0, 0)/e^ ^ TraceG(o,o,2)d2gTrace(nj._G-)x 



(4.5) det 



u'i-x)\Vx ••• 7^- 7^+ o ••• o 



evaluated at xi = 0, , A) = (0,0), where V'j := bii4>j, ( ) is the part of the 

\\U [-X)\VxJ 

first column of ¥~a involving only slow modes 

(recall (|2.11|) ). and ^•^-^ part involving only slow modes of the block of + ¥~a 

corresponding to the slow block ^•j^-^ °^ ^ • Referring to (j4.2p . we see that Vx has 

a limit V as X ^ +00 so long as S := lim2^-oo(f^VI^'l)(^) exists (the first column 
of {—X)e~'^ being then approximately \U'{—X)\S as X ^ +cc, so that the slow 

component of the first column of F~a is approximately |?7'(— X)|y for a fixed V determined 
by S). 

By a block determinant expansion, we have, therefore, 

Dx{0, 0) /e-^o""^ TraceG(0,0,2)dZgTrace{nj._G-)X ^ 

a\U'{-X)\ det(V2, . . . , V'n+i) det(7^-, 7e+, Vx)L,=o, (€,A)=(o,o) 
^ a\U'i-X)\ det(V^2, . . . ,0n+l)det(7^-,7^+,y)|^.^^o_(^-;^)^(o_o), 

where a = ±1 depending on dimensions n, r and det(02) • • • ,'4'n+i) assuming stability 
of the limiting shock U (else P would vanish to second instead of first order at the origin 
|Z3| ). Normalizing 

sgn(Tdet(V'2, • • • , V'n+i) = +1, 



we obtain 



lim sgn lim Dx{0, X) = 6 := det{R' , R+,V) 
A^o+ real 



as claimed. 

Reviewing the computation in [Z3] of dx'D{0, 0), we find that this is the same normaliza- 
tion of ipj columns leading to the assumed normalization (|1.14|) , whence (|1.15|) is necessary 
for stability by the discussion above the statement of the theorem. 
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Finally, for a Lax 1-shock, all eigenvalues of are positive, hence the TZ'^ block is 
empty in the computation above and the boundary conditions must be of inflow type. 
If also, the boundary conditions are Dirichlet type, then by (|2.3p . the first column of 
E0(0, 0)6"'^"*'^"'^'-^ (since all of E^) is of form 



1 




where E-,^ and E^ defined as the first column of E := T {—X)e '''"E and the 
slow part of F~ := T~ {—X)e~'^ ^¥~, are asymptotically of form (since T~e^ ^^E^ = 
(6iiJ7'(xi), 0) and lim^^-ooT^ = Id) 

\U'{-X)\ (^'f^ and \U\-X)\ {^''^'^f^') 

as X — > +00, where TZ^ lies in the unstable subspace of . Equating, we find that 
a- = and h^-^S = h^^{A];)-^n:^ . Recalling that {Af.Af^U' = by the linearized 
boundary-layer ODE, so that 

({A-,f\{A-,)'^)s = Q, 

we find by inspection that IZi = A^S, or, in the notation of the general case, 
(4.6) V = ^^5, 

yielding the result. □ 

Proof of Theorem \1.13[ Immediate, by Lemma 13.21 and Theorem 11.111 □ 

Proof of Theorem \1.15i Set r] := \U'{—X)\, and consider < p < Ce^^^ for some fixed 
C > and 9 > sufficiently small, in particular, small enough that 



(4.7) \U'{-X)\«Ce 



-ex 



setting 

(I, A,r/) =: (pIcpAo, Wo) 

with 1(^0) ''^O) %)| = 1- Note, as X — > +oo, that the set of possible values of (Coj^o^Vo) 
expands to the full positive half-sphere ^Xq, rjo > 0. Restrict now to a compact subset of 
the positive half-sphere, recalling (see Proposition 12. 2[) that the Evans functions Dx and 
V, and their component columns, are C°° in p, ^q, Xq. 
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Within the specified parameter-regime, both slow and fast modes of (j2.ip at A) = 
plCoi-^o) are well-approximated on xi G [— X, 0] by their limiting values as p ^ 0, de- 
scribed in Proposition I2.4[ Mimicking the one-dimensional computations (j4.5|) . (|4.6p . we 
may rewrite 



Dx{pio,p\o)/e^o TraccG(pCo,pAo,.)d, 



Iz Trace(n„_G-)X 



as the sum of 



and 



pdet 



^^^*V%^i'5 ••• 7^- 7^+ ••• o 

+ o{pr]o) 



dpT-E^ - dpT+El ■■■ -R- TZ+ • • • 

+ o(p|(lo,Ao)|), 



both evaluated at xi = 0, p = 0, and {(,o,Xo)- We omit the details of this straightforward 
but tedious computation. 

A standard computation |ZS[ IZ3[ IGMWZ4] using the variational equations of (|2.ip with 
respect to p yields 

d 

dpT-E^ - dpT+Et = Xo[U] + Y^iil[F,{U)], 

i=2 

whence, normalizing as usual so that 

sgncjdet(V'2, • • • jV'n+i) = +1, 
we obtain by block determinant expansion 

p-^Dx{pio,pXo)/e^O~'' TraceG(p|o,pAo,z)dZgTrace(n^_G-)X ^ 

A(e~o,Ao,r/o) + o(l), 

where A is defined as in (jl.2ip and o(l) is with respect to Xq, and tjq for each fixed 
X and — > uniformly as X ^ 0. 

By an application of the Implicit Function Theorem, it follows that existence of a root 
(^Q, Aq, ryg) of A on JRAq > at which dx^A ^ implies existence of a nearby root (fg, Aq, ?7q), 
3f?Aj > 0, p > 0, of 



2 Trace(n„_G-)X 



hence of Dx, for X sufficiently large, or instability of U'^ . Thus, nonvanishing of A on the 
strictly positive half-sphere is necessary for stability as X — > 0. □ 
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Proof of Theorem \l.l(A The estimate (|4.8p in fact holds for all r/ := \U'{—X)\ and < p < 
Ce~^^, with the o(l) term uniformly decaying and uniformly C*^ as X ^ oo (however, 
not uniformly see |GMWZ4l IHMWZSIIUMWZG] ^. It follows therefore, that nonvanish- 
ing of A on the (closed) nonnegative half-sphere, implying a lower bound on |A|, implies 
nonvanishing of Dx on the parameter range < p < Ce-^^, for X sufficiently large. 
If Ce~^^ < p « 1, on the other hand, a much cruder estimate yields 

p-l Dx{pio,p\o) /e^o^ TraceG(pfo,pAo,2;)a!2gTrace(ng_G-)X ^ 

A(|o, Ao) + 0(1) + 0{\U'{-X)\/Ce-''') = 
A(eo,Ao) + o(l), 

by ()4.7p . again with o(l) uniformly decaying as X — > +00. This implies nonvanishing of 
Dx on the parameter range 1 >> p > Ce~^^ , for X sufficiently large. 

For p bounded from below, on the other hand, we have by the basic convergence result of 
Theorem 11.81 that Dx is nonvanishing if V and D_ are nonvanishing, i.e., if U and U = U- 
are stable. This completes the proof of the theorem. □ 

5 Application to gas dynamics 

We now apply our results to the fundamental example of compressible gas dynamics, re- 
stricting without loss of generality (by rotational invariance of the equations) to dimension 
d = 2. Consider the compressible Navier-Stokes equations 

(5.1a) pt + {pu)x + {pv)y = 0, 



(5.1b) {pu)t + {pU^ + p)x + ipuv)y = (2/U + r])Uxx + PUyy + {p + r])Vxy, 



(5.1c) 



{pv)t + {puv)x + {pv^ + p)y = pvxx + (2/i + r])vyy + {p + ■q)uyx , 



(5.1d) [pE)t + {puE + up)x + {pvE + vp)y 

= {kTx + {2p + r])uux + pv{vx + Uy) + rjuv^ 

+ ( KTy + {2p + ri)vVy + PU{VX + Uy) + TJVUx 



on the half-plane x G M^, y € M, where p is density, u and v are velocities in x and y 
directions, p is pressure, 



(5.2) ^ = ^ + T+2 
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is total energy density, e and T are internal energy density and temperature, and constants 
/i > \r]\ > and k > are coefficients of first ("dynamic") and second viscosity and heat 
conductivity. 

We assume ideal ("7-law") gas equations of state 

(5.3) po{p,T)=Tpe, eo{p,T) = c^T, 

where c^, > is the specific heat at constant volume, F := 7 — 1 > 0, and 7 > 1 is the 
adiabatic index of the gas; equivalently, 

(5.4) p(7;,5) = ae^/'="/9^, 

where S is thermodynamical entropy \Ba\ ISmj . In the notation of (jl.2p , we have 

(5.5) U = {p,pu,pv,pE) =: {p,mi,m2,8) 
and 

(5.6) = {pu, pu^ + p, pv"^ , puE) = {mi,m\/ p + p,mim2/ p,mi£/ p). 

Remark 5.1. In the thermodynamical rarified gas approximation, 

2n + 3 97-5 2 

(5.7) ^ = ^' -/^=^ ^ = -3^ 

for V := «;/c„, where n is the number of constituent atoms of gas molecules (here assumed 
to have "tree" structure) |Baj . with 7 = 5/3 and 7 = 7/5 for the main applications of 
monatomic and diatomic gas. In particular, 

(5.8) 1< 7 < 2 and — - — > 1 

2p + r] 

for common gases, a conclusion that is born out by experiment. See Appendices A and B 
of [HLyZlj for further discussion. 

5.1 Viscous Shock Profiles 

From (15. ip . setting time-derivatives to zero, integrating in x, and rearranging, we obtain 
after a brief calculation the standing-shock ODE 



v! = (2/i + r?) ^ ^m(n — ti_) -|- T{pe — />_e_)^ , 



(5.9) .2 

-/ -i( I" \ m[u — u^) 



e' = V " ( m{e — e_) ^ \- {u — U-)T p^T^j , 



where m := pu = constant and v = constant. 
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Using various scale- invariances of system (jS.ip . we may take without loss of generality 
m = p- = = 1, = = 0, yielding 

^' = -^((n-l)+r(4-e- 

{u-iy 



e' = u- 



M(e-e„)-^^^-^ + (n-l)re- 



with V = 0, with endstates 

u+a(u+ - 1) {u+ - l)(r + 2) 

^+ = r(r + 2-a) ' r(r + 2-a) ' ^+ = '/"^ 

a := ^^^Z^""*" 1 parametrized by the single quantity 

T 

(5.12) 1 > n+ > := 



r + 2 



In the strong-shock limit — > u^,, e_ ^ 0, with all other quantities remaining in physical 
range; for details of these computations, see |HLyZ2| , Sections 3-5. 
Linearizing ()5.10p about (n_,e_) = (l,e_), we obtain 

(-) (:)'--(:)■ --(tDC-"' 

determining the asymptotic behavior of (n, e)(z) as z ^ —oo. One may check for all 
1 > u+ > that M_ has two positive distinct real eigenvalues 0<a;_<l/i^< a;+, 



l-Fe- _ 1 \ I // 1-re- _ 1\2 4re- 



with associated eigenvectors sj = (— 1, — ^^(^^^1"^/^) merging in the special limiting case 
u+ ^ ti*/e_ ^ 0, 2/x -|- 7/ = to a pair of real semisimple eigenvalues. 

That is, for a Lax 1-shock, [/_ is a repellor for the standing-wave ODE, and f7+ a saddle, 
in agreement with the abstract conclusions of |MaZ3j for extreme shocks of general systems 
and of [Gi] for shock profiles of gas dynamics with general equation of state. In particular, 
note that 

detM_ = {u{2n + r]))-\l-r{l + T)e^) > 0, 
with (1 — r(l -|- r)e_) approaching 1 in the strong shock limit u — > ti*/e_ 0, and 
^ ~ (r+2)^^ ^ (r+2)^ weak shock limit ?x+ ^ 10 By reality and simplicity of the 

eigenvalues coj, we have that limits 



(5.14) s:= lim {u' ,e')/\{u' ,e')\) 

- ■ oo 



^ This repairs an error of iSZj, in which U- was mistakenly computed to be a saddle, leading to an 
incorrect value of S. 
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and 



(5.15) 



S := lim {U'/\U'\) 



I -1 0\ 

dU _ 

d{u,ey^-^~ 

Vl/2-e- 1/ 

exist, with s generically lying parallel to the slow mode s_, or 



S 



(5.16) 



-1 

\ 



(5.17) 



/ -1 o\ 



Vl/2-e_ ly 

/ 1 




V~v{u)--l/v) _ 1 



Finally, from ()5.5p - ()5.6p . we obtain after a brief calculation 

/ 1 0\ 

Pp-l+Pe/2 2-pe p, 

10 

V -1/2 1/2 1/ 

/ 1 0\ 

re_-i + r/2 2-r/2 o r 

10 

V -1/2 1/2 1/ 



■.= d{Fi/dU){U-] 



from which we compute 



(5.18) 



A];s 



re.-i + r/2 + re_ ^;g--;y 



^_ r-u{uj.-i/u) _ ^ 

V u{uj- — l/u) 



(5.19) 



5.1.1 The strong shock limit 

2/7T»? ^ ^^-^^P converges to 5 = (1,0,0,-1/2)-^ in the strong shock limit e_ 0. 
For jjTPn — however, ()5.16p becomes singular in the limit as e_ — > 0, for which also 
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LO- 1/u. To evaluate this limit, it is easier to return to (j5.13p and compute directly with 
e_ = 0, to obtain s — > (—1, 1 — 0)^, yielding the general formula 



(5.20) S (1,0,0,1/2 -min{l,(/)})'^ as tt+ ^ n_, 

Noting that 



2fi + ri 



(5.21) 



we thus have 








1 





0\ 




-i + r/2 


2 - r/2 





r 










1 





v 


-1/2 


1/2 





1/ 



(5.22) 



( \ 

r max{0, 1 — </)} — 1 



V -min{l,0} / 



and 
(5.23) 

completing our asymptotic analysis. 



[f/]-(pO,0,^^" 



5.2 One-dimensional stability 

Following the treatment in [Set ISZ| , we note for Lax 1-shocks that 
(5.24) det(7e+,/) =^+-/, 

for any vector / G C", where is the unique stable left eigenvector of A+ and • denotes 
complex inner product. In one dimension, is just the stable left eigenvector of which 
may be computed to be 



(5.25) 

where 
(5.26) 



Pp + CU + 



Pe{v?/2-e) PeU 



P 



P 



c,0,^) {U, 



c+n_i_ + rti^/2, — rn+ — c+, 0, r ) , 



c : = 



denotes sound speed. This computation is most easily accomplished by working in the 
more convenient nonconservative coordinates {p,u,v,e), which are related to conservative 
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variables {p, pu, pv, p{e + + /2)) by a readily computed lower triangular change of 
coordinates; see or Appendix [Xl 
Combining all facts, we have 

5 = l+-A^S 



(5.27) - V - - V'- - 1 + + re- _ r/,) 

/ r - v{u^ - i/u) 

It is readily verified on the other hand that 

(5.28) (5 > 0; 

see Section [5.2. II just below. The one-dimensional stability condition (|1.17p thus reduces in 
this case to 

(5.29) sgn^ > 0, 

a condition that can be readily checked numerically using (|5.27p . 

5.2.1 The strong shock limit 

In the strong shock limit — s- n^,, we have e_ — > 0, a — > +00, and 

e+ ^n,(l-u,)/r = 2/(r + 2)2, 



so that c+ y2f(r + l)/(r + 2) and 



rV2r(r + 1) _ rV2r(r + 1) y 

(r + 2)2 +2(r + 2)2' r + 2 r + 2 J ' 



hence, by 



5 = ^+■ A^S 



r max{0, 1 —(/)} — 1 1 — r min{l, 0}. 



Meanwhile, 



r+2 r+2 



from which we may conclude by homotopy/nonvanishing of 5 that 5 > for all 1 > u+ > n^=, 
verifying dOHjl - dOgl ). 
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The case (/) > 1. For </>>!, (|5.29p becomes 



r./2r(r + 1) n 

— — ^ — ■ — - H 1 - r > 

r + 2 ^r + 2; 

or 

V2r(r + 1) > 2, 

which evidently fails for T in the kinetic range < F < 1 (indeed, for all T outside (1,2)). 
Thus, we may conclude instability in the strong shock limit in this range. 
The case < 1. For (j) < 1, (|5.29p becomes 



Fv^2F(F + l) 



F+2 F+2 
or 



F(l - 0) - 1 - F(/)> 



V2F(F + 1) + F 1 - F(l - (/.) -(F + 2)0>O. 



Defining a := y^2F(F + 1) + F, we may rewrite this as 

C7(l -F) > {T-Ta + 2)<j), 

or, assuming F(l — o") + 2 > 0, as holds for example on the kinetic range < F < 1, or 
1 < 7 < 2 (on which c7<2 + F, so2 + F> Fa), as 

(5.31) 0< 



F(l-a) + 2' 

which is satisfied for small enough, but for (j) = 1, hence for (f) < 1 large enough, is not 
satisfied, by the analysis of case = 1 above. 

Common gases and the kinetic approximation. Recall that for common gases, (p 
is less than one. For gases obeying the kinetic approximation (|5.7p - (|5.8p . 

(5.32) 4 



27F + 12 ' 



so that (j) < 1 for F > 4/27 ~ .148, in particular for n-atomic gases with n < 5. Thus, it is 
the case (p < 1 that is relevant to typical applications. Substituting (]5.32p into (15.3ip and 
noting that y^2F(F + 1) < F + 1 for < F < 1 yields the necessary condition 

(5.33) 16(F + 2) < (2F + 1)(1 + 15F), 

or < (F — 1)(30F + 31), which is violated for the entire kinetic range < F < 1. 

Conclusions By Remarks 11.121 2 and ll.l8"l 2. boundary layers are both one- and multi- 
dimensionally stable in the standing shock limit for limiting shocks of sufficiently small 
amplitude, i.e., 1 — tt+ sufficiently small. By the calculations above, however, for typical gas 
laws, they are not even one-dimensionally stable in the strong shock limit for limiting shocks 
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of sufficiently large amplitude, i.e., u+ — u^ sufficiently small, even though the corresponding 
shock is perfectly stable [HLyZl HLyZ2| . 



Thus, we have the striking conclusion that for (all!) typically physically occurring gases 
under inflow Dirichlet boundary conditions, there is a transition from stability to instability 
of boundary layers in the standing shock limit as the amplitude of the limiting shock increases 
from zero to its maximum value. 

5.3 Multi-dimensional stability 

The computation of ^+(^,A) in multi-dimensions may be found, for example, in Appendix 
C, |Z3pl . where it is computed as 

(5.34) l+{i,\) = [e ; +~r' / ---vu,—^^=,r] 

2 _ «2 P Jc2_P2 Jc2_a2 



where O ■= Pp — ^ = 2Te, r] := ^ =T, and c is sound speed (|5.26p . or 



(5.35) 4(1. A) = (2re, - ^^±S±^ + ^< - Tu,, T 

where 



(5.36) P :-- 



— 



Together with our computation of A^S in (|5.18p . this determines • A\S. Meanwhile, 
A := ■ (A[f7] + z^[F2(f7)]) is computed for the same choice of 1+ in Appendix C, |Z3j 
(equation displayed below C.36), thus determining ?7(^, A) = • A\S/ /S.{^,X). 

With Proposition 11.191 this gives a straightforward means of numerical determination 
of multidimensional instability, by plotting the image of 17(1, ir) as r ranges over the real 
axis and checking whether or not this curve strikes the nonnegative real axis; however, we 
shall not carry this out here. 

The numerical determination of one- and multi-dimensional stability transitions for ideal 
and other gas laws would be interesting problems for further investigation. A further very 
interesting open open problem is to determine analytically the stability transitions as was 
done for the inviscid shock problem (involving only A) in |Erl [M] ; see Appendix C, 



A Computation of in one dimension 

In this appendix, we carry out for completeness the computation of for the one-dimensional 
Navier-Stokes equations, verifying (I5.25p . In variables (p, u, u, e), the quasilinear hyperbolic 

Contributed by K. Jenssen and G. Lyng 
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part of the equations becomes 

pt + q- Vp + pdivg = 0, 
(A.l) qt + q-Vq + p-^PpVp + p-^pe^e = ^, 

et + q ■ Ve + p^^pdivu = 0, 

where q = (u, v) denotes velocity, or, in one dimension, 

Vt + {uld + M)F,, = 0, 

where V = {p, u, v, e) and 



M :-- 



/ p \ 

p'^Pp p^^p, 



V p-^p y 



from which we may conclude that Ai = S{uld + M)S ^ for 

/ 1 



d{p, pu, pv, p{e + + v'^/2)) 
d{p,u,v,e) 

( 



0\ 

u p 

V /9 

\e + + 't;^/2 pu pv p) 



V 



1 








o\ 


p 

— V 

p 


1 
p 




1 

p 






e+u^/2+v'^/2 


— M 


—11 




P 


p 


p 



and thus = ^llj.S'^^ for ^+ defined as the left eigenvector of M associated with the 
eigenvalue of smallest real part, * denoting adjoint, or congugate transpose, all quantities 
to be evaluated at {p,u,v,e) = (l/n+, 0, e+). 

By inspection, = {pp, —pc, 0,Pe) for v = 0, where sound speed c is defined as in (|5.26p . 
whence 



Pp + CU + 



PeW'l^ - e) PeU 



c,0,^)(C/4 

p 



as claimed. 
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